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Recent astronomical observations indicate that the Universe is presently almost flat and undergo- 

\ j^ , ing a period of accelerated expansion. Basing on Einstein's general relativity all these observations 

(— ^ ' can be explained by the hypothesis of a dark energy component in addition to cold dark matter 

("^ ' (CDM). Because the nature of this dark energy is unknown, it was proposed some alternative sce- 

^SJ . nario to explain the current accelerating Universe. The key point of this scenario is to modify the 

I ■ standard FRW equation instead of mysterious dark energy component. The standard approach to 

'~{ ' constrain model parameters, based on the likelihood method, gives a best-fit model and confidence 

^~~i , ranges for those parameters. We always arbitrary choose the set of parameters which define a model 

(/-v ■ which we compare with observational data. Because in the generic case, the introducing of new 

fv^ ' parameters improves a fit to the data set, there appears the problem of elimination of model pa- 

\ rameters which can play an insuflicient role. The Bayesian information criteria of model selection 

^\l . (the AIC and BIC) are dedicated to promotion a set of parameters which should be incorporated 

^ ' to the model. We divide class of all accelerating cosmological models into two groups according to 

f~~« ' the two types of explanation acceleration of the Universe. Then the Bayesian framework of model 

fSJ , selection is used to determine the set of parameters which gives preferred fit to the SNIa data. We 

CO ' find a few of flat cosmological models which can be recommend by both the Bayes factor and Akaike 

^^ ' information criterion. 

o: 
^. 

O ■ I. INTRODUCTION 

Oh' Recent measurements of distant supernovae (SNIa) [ij, |3| as well as current measurements of cosmic microwave 

Q , background anisotropics Q favor a spatially flat Universe filled by cold dark mater (CDM) [J, Q and a dark energy 
$H ' component of unknown origin [l|, Q] , which can be modelled as a perfect fluid with energy density px and pressure px 
c/3 , violating the strong energy condition px + 3px > 0. The combination of CMB and SNIa data with other orthogonal 
Cu ' measurements, as the HST determination of the Hubble parameter, constrain the Universe to be almost flat even if 
^ \ we consider the time variation of dark energy equation of state |6j . All these models explaining current acceleration 
of the Universe in terms of smoothly distributed and slowly varying 'dark energy' are formulated in the context of the 
standard cosmological picture based on general relativity theory. Possible types of explanation include a cosmological 
constant K ^, an evolving scalar field (quintessence) [8| , the phantom energy |9|, Il0| in which a weak energy condition 
is violated, models filled with Chaplygin gas [ll|, models with dynamical coefficient equation of state wx = Px/ Px 
(decaying vacuum energy density), usually parameterized by the scale factor a or redshift z, where (1 + z = a~^), fluid 
which describe quantum effects coming from massless scalar field (the Casimir effect) [l3;ll3;ll3; the noninteracting 
Chaplygin gas and baryonic matter [isL Ha . [l7| . In Table |l] we collected the ten candidates for dark energy description 
together with their Friedmann first integral 

3ij2=Peff-^, (1) 

where H ~ {dhia)/{dt) is the Hubble function, and fc = 0, ±1 is the curvature index). For all these approaches some 
hypothetical dark energy component is postulated which satisfies the conservation condition 

Pi = -3H{pi+p^), i = m,X (2) 

for both standard dust matter as well as dark energy X separately. 

On the other hand, alternative ideas to the dark energy idea have been recently offered. Freese and Lewis [ig 
proposed so called the Cardassian model, in which the Universe is flat, matter dominated and accelerating not due 
to dark energy but as a consequence of modiflcation of the Friedmann first integral as follows 

3H^^p + Bp", (3) 



^ 



TABLE I: The Hubble function for 10 prototypes of cosmological models explaining the present acceleration of the Universe in terms of dark energy 
case model _ff^ (2;) relation independent model parameters 

1 ACDM model H'^ = Hi{nm,o{l + zf + {1 - ^^,0)} -ffo, fim.o 

2 non-flat FRW model with A H'^ = Ho{n^,o{'i- + zf + ^kA'^ + ^f + ^J^,o} Ho, Qm,(), ^a,o 

3 FRW model with 2D topo- H'^ = Hi{nn,,o{l + zf + 0^,0(1 + zf + aopo(l + z)} Ho, fim.o, fife.o 
logical defects px ~ ^%Px 

4 FRW model with phantom H'^ = ffo{^m,o(l + zf + Q,k,a(l + z)^ + ^phs)(l + z)'^} Ho, ^mj,, fifc,o 
dark energy px = -^px 

5 FRW model with phantom H'^ = Hl{n^,o{l + zf + t^fc,o(l + zf + !^ph,o(l + «)3(i+"'^)} Ho, ^^,0, Ofe.o, wx 
dark energy px — wxpx, 

Wx < — 1 fitted 

6 FRW model with Chaplygin H^ = H^ [n^A^ + ^f + "fc.o(l + zf + f^ch,o[^s + (1 - As){l + zf]H Ho, n,n.o, O^.o, As 
gaspx = -^, ^>0 

7 FRW model with general- H^ = H^ S^n^A^ + zf + ^kA^ + ^f + ^Ch,o[As + {1 ^ As){l + zf'-'+"^]Th^ Ho, ^,,.,0, ^k,o, As, a 
ized Chaplygin gas px = 

— 4-, A > 0, a = const 

8a FRW models with dynami- H'^ = H^ l^mA^ + zf + ^^,0(1 + zf + t7x,o(l + z)3(»o-»i+i) exp[3«;i2]| Ho, f^m.o, f^x.o, wo, wi 

cal E.Q.S. parameterized by 
z Px = (wo +wiz)px 

8b FRW models with dynami- H^ = i^o^ [n^A^ + zf + f^fc.o(l + zf + !^x,o(l + 2)3(-o+-i + i) exp[-2^]} /fo, an,o, ^x.o, wo, wi 

cal E.Q.S parameterized by 
scale factor a px ~ (wo + 
Wi(l - a))px 

9 FRW model with quantum H^ = _H'o{^m,o(l + zf + i^k,(){l + zf + Qa,o - ricass,o(l + zf} Ho, rJm,(), nA,o, ^Cass.o 
effects origin from massless 

scalar field at low tempera- 
ture (Casimir effect) px = 

^3-, PX,0 > U 

10 flat FRW model with Chap- H^ = Hi Ul - n^AlAs + (l - As)il + zf^i-"'^]^^ +nA'^ + zf] Ho, As 
lygin gas and baryonic mat- 
ter peff = - Sap^'H; 

f2m,o = 0.05, m = —1.5 are 
fixed 



where B is a constant and the energy density contain only dust matter and radiation (SttG = c — 1). The additional 
second term in relation (O rnay arise from 'new physics'. It dominates at the late epoch and drives the present 
acceleration of the Universe. Because terms of type p" may arise as a consequence of living on (3 + 1) brane in a 
high-dimensional bulk space, the origin of acceleration lies rather in modification of the FRW equation at low energy 
scales than due to a dark energy component. In brane world scenarios, the observer is embedded on the brane in 
a larger space on which gravity can propagate. The idea is that an observer measures only 4-dimcnsional gravity 
up to some corrections that given the weakness of gravity, can in gen eral be made small enough not to conflict with 
observations without tweaking with a model parameter too much [l9j . 

In Table E] it is collected ten representative models offering explanation of current acceleration of the Universe 
in an alternative way to a dark energy. Apart from the Cardassian model there are different brane world scenarios 
which were originally proposed by Dvali, Gabadadze and Porrati (DGP) [23|, Deffayet, Dvali and Gabadadze (DDG) 
|21| , Randall and Sundrum 22] , Shtanov brane models [23| • We also consider models of 'nonlinear gravity' based on 
modified Lagrangian density which is proportional to i?", where R is the Ricci scalar and n is constant |2J| (n = 1 for 
standard general relativity) , models based on non-Riemannian gravity |23 , bouncing models arising in the context of 
loop quantum gravity and models with energy transfer between the dark matter and dark energy sectors (A decaying 
vacuum and phantom dark energy). 

The main goal of this paper is to make the ranking of accelerating models using the Bayesian framework |26t |27| . 
The effectiveness of application of information criteria of model selection in the cosmological context was recently 
demonstrated by many authors 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44). In the paper we obtain 
these accelerating models which are the best ones from the set under consideration in explaining the SNIa data. We 
calculate the Bayes factor for all models with different numbers of parameters, which differentiate between them '455 • 
The Bayes factor measures the change in relative probabilities of any two models in light of observational data (SNIa 
data) when we update the prior relative model probabilities to the posterior relative model probabilities. 

In observational cosmology many theoretically allowed models with different prediction of the past (big bang versus 
bounce) and the future (big rip versus de Sitter attractor) becomes in good agreements with the observational data 
pq. Therefore we propose to use the Bayesian factor to differentiate among all these models and make some ranking. 

II. IDEA OF MODEL SELECTION 

In the development of cosmology the basic role played an idea of cosmological models together with an idea of 
astronomical tests l47|. The idea of cosmological tests make cosmological models parts of astronomy which offers 
possibility of observationally determining the set of realistic parameters, that can characterize viable models. While 
we can perform estimation of model parameters using a standard minimization procedure based on the likelihood 
method, many different scenarios are still in a good agreement with observational data of SNIa. This problem 
which appears in observational cosmology is called the degeneracy problem. To solve this problem it is required to 
differentiate between different dark energy models. We propose to use the Akaike information criterion (AIC) and 
BIG quantity (as an approximation to the logarithm of the marginal likelihood). 

For the model selection framework it is required to have data and a set of models and then we can make the model 
based statistical inference. The model selection should be based on a well-justified single (even naive) model or, at 
least, a simple model which suffices for making inferences from the data. In our case the ACDM model plays just 
the role of such a model. The model selection should be viewed as a way to obtain model weights, not just a way to 
select only one model (and then ignore that selection occurred). Moreover in the notion of true models do not believe 
information theories because the model by definition is only approximations to unknown physical reality: there is no 
true model of the Universe that perfectly reflect large structure of space-time, but some of them are useful. 

In this paper the Bayes factor and AIC quantity are used to compare models gathered in Tableland Table Ull 

The AIC is defined in the following way |48l | 

AIC = -21n£ + 2d, (4) 

where C is the maximum likelihood and d is a number of model parameters. The best model with a parameter set 
providing the preferred fit to the data is that which minimizes the AIC. While there are justification of the AIC in 
information theory and also rigorous statistical foundation for the AIC, it can be also justified as Bayesian using a 
'savvy' prior on models that is a function of a sample size and a number of model parameters. For the AIC we can 
define AAICi as the difference of the AIC for model i and the AIC value for the reference model: A^ = AICi — AICi. 
The Ai are easy to interpret and allow a quick 'strength of evidence' comparison and a ranking of candidates for dark 
energy description. The models with < A^ < 2 have substantial support (evidence), those where 4 < A^ < 7 have 
considerably less support, while models having A^ > 10 have essentially no support with respect to model 1. 



TABLE II: The Hubble function for 10 cosmological models beyond the standard general relativity 



case model 



H (z) relation 



independent model parameters 



1 Cardassian type of Fried- 
mann equation, Q,r,{) ~ 10^'* 
is fixed 

2 Dvali-Gabadadze-Porrati 
brane models (DGP) 

3 Deffayet-Dvali-Gabadadze 
brane models with A (DDG) 

4 Randall- Sundrum brane 
models with dark radiation 
and A = 

5 Randall- Sundrum brane 
models with dark radiation 
and A (RSB) 

6a Shtanov brane models 
(Branel) 

6b. Shtanov brane models 
(Brane2) 

7 modified afiine gravity 
(MAG) model 

8 FRW models of nonlin- 
ear gravity with Lagrangian 
density proportional to Ricci 
scalar R (NG) 

9 bouncing models with A 
(BACDM) 

10a models with energy transfer 
(dark matter <^ vacuum en- 
ergy pi sector) 

10b models with energy trans- 
fer (dark matter ^-> phantom 
dark energy sector) 



H^ = Hi I Q.k,o{l + zf + fin.,o(l + zf 



TTI + (1 + ^) 



S2m,0 



1+z • n„.o 
. "^ n' 



H^ = H, 



(?{ [7^^,0(1 + 2):^ + 0^e,() + V'fW] + f^fc,0(l + zf} 



H'^ = Hi 



^ + Jf^m,0(l + 2)^ + fiA.O + J-hr? 



2roHo 



H^ = HiiQmA^ + zf + «fc,0(l + zf + ftdrfii^ + zf + !^A,o(l + zf} 

H"" = Hi{n„,,o{l + zf + fifc,o(l + zf + fid,,o(l + zf + !^A,o(l + zf + fiA,o} 

H^ = Hi 1^^,0(1 + zf + n^fi + 2nifi - 2ynZoV'^m,o(l + zf + n<,,(, + «!,(, + QAbfi] 

H^ = Hi |t7m,o(i + zf + n^,o + 2nifi + 2^nf^^nn,,o{i + 2)3 + ^,,,0 + 0.1,0 + fiAb.o} 

H^ = Hi{nm,o{l + zf + fifc,(,(l + zf + fiv,,o(l + zf + ^A.o} 



H^ = i/o'{«m.o(l + Zf^ + n,,o(l + zfi^^^}Qnonl,o{^ + z) 



3(l-'») 



{n-3)2 



H^ = Hiinn.A'^ + zf + !^fc,o(l + zf - n„,o(l + zf + Oa,,)} 

H^ = Hii^mA^ + zf + ant(l + zf + Oa.o} 



H^ = H 



i {t7n,.o(l + zf + ant(l + zf + nph(l + 2)3(1 + "'^)} 



Ho, fife.o, n 



Ho, fim,0, firc,0 



Ho, To, Oa.O 



Ho, Qmfi, fidr.O, f^A.O 



Ho, fim.O, fidr.O, ^A.O, ^A.O 



Ho, ^1X1,0, fiiT.o, fi;,o 



Ho, fim.O, Oa.O, Ol.O 



Ho, f2m,0, ^1^,0, f^A.O 



Ho, ^1X1,0, Hnonl.O, « 



-ffo, fim,0, fin,0, Oa.O 



^0, fim.O, Oa.O, 



Ho, Slm,0, fiph.O, Wx 



l(^ 



In the Bayesian framework a best model (from the model set {Mi}, i — I, . . . , K) is that which has the largest value 
of probability in the light of data (so called a posteriori probability) . We can define the posterior odds for models 
Mi and Mj , which (in the case when no model is favored a priori) is reduced to the marginal likelihood (E) ratio (so 
called the Bayes factor - Bij) 

^ _ jL{e\D,M,)P{e\M,)de ^ e, 

'' jL{f]\D,M,)P{f]\M,)dv E,' ^> 

where ^ is a parameter vector, which defines model i, L{d\D,Mi) is likelihood under model i, P{6\Mi) is the prior 
probability for 9 under a model i. 

Schwarz |49j showed that for iid observations coming from the linear exponential family distribution the asymptotic 
approximation {N — > oo) to the logarithm of the marginal likelihood is given by 

ln£; = ln/:--ln7V + 0(l). (6) 

According to this result he introduced a criterion for the model selection: the best model is that which minimized 
the BIC, defined as 

BIG = -21n/: + dlniV. (7) 

This criterion can be derived in in such a way that it is not required to assume any specific form for the likelihood 
function but it is only necessary that the likelihood function satisfies some non-restrictive regularity conditions. 
Moreover the data do not need to be independent and identically distributed t>G\. 

To compare models Mi and Mj one can compute 2\nBij — —(BIC; — BICj) = — ABICy which can be interpret 
as 'strength of evidence' against j model: < 2ln Bij < 2-not worth more than a bare mention, 2 < 2\n Bij < 6 - 
positive, 6 < 2\n Bij < 10 - strong, and 2ln Bij > 10 - very strong. 

It is useful to choose one model from our model set (a reference model-s) and compare the rest models with this 
one model, situation in which 21ni?si > indicates evidence against model i with respect to the reference model, 
whereas 2\TiBsi < denotes evidence in favor of model i. 

We can compute posterior probability for model i in the following way 

PiM,\D) = /'^ , (8) 

where Bis = exp[iAi3/Csi]. Then one can see how prior believe about model probability P{Mi) ~ j^ change 
after inclusion data to analysis. This is the probability for model i being the best model from set of models under 
consideration. 

There are many simulation studies in the statistical literature on either the AIC or BIC alone, or often comparing 
them and making recommendation on which one to use. It should be pointed out that both of them are an asymptotic 
approximation. This assumption is satisfied when sample size used in analysis is large enough, large with respect to 
the number of unknown model parameters. 

Note that the assumptions of using model selection methods are satisfied, namely 

1. there is model-based inference from SNIa data (the Imiiinosity distance observable); 

2. there is a set of models and no certainty about which model should be used in explanation of present acceleration; 

3. a data-based choice is made among these competing models (see Tabled Table Hljl . 

In this paper the restricted 'Gold' sample of A^ = 157 SNIa [l| is used. It is assumed that the supernovae 
measurements come with uncorrelated Gaussian errors. In this case the likelihood function has the following form 

-1 / -^ /..thcor ,,obs\2\ 

i^ocexp 1 V (/-^ -/^ ) , (9) 



where (Ji is known, fi°^^ = rrii — M (mi-the apparent magnitude, M-the absolute magnitude of SNIa), ^f^^°'^ — 
5 logio ^Li + A^, M = —b logio -^0 + 25 and D^i = HodLi, where dLi is luminosity distance 

dL^^(l + Z,)-^^l=T(HoJ\n^\ [^ -^), (10) 

^0 V?M V ^ Jo H{z')J 



TABLE III: Values of the AIC, BIC, AAICis and 2 In Bsi (where s is the index of the reference model) for the flat models from 
Table n 



case 


AIC 


BIC 


AAICii 


21nBH 


AAICio 


2 In Bio, 


1 


179.9 


186.0 


0.0 


0.0 


2.0 


2.1 


2 


179.9 


189.0 


0.0 


3.0 


2.0 


5.1 


3 


183.2 


189.4 


3.3 


3.4 


5.3 


5.5 


4 


178.0 


184.1 


-1.9 


-1.9 


0.1 


0.2 


5 


178.5 


187.7 


-1.4 


1.7 


0.6 


3.8 


6 


179.7 


188.9 


-0.2 


2.9 


1.8 


5.0 


7 


181.7 


193.9 


1.8 


7.9 


3.8 


10.0 


8a 


180.5 


192.7 


0.6 


6.7 


2.6 


8.8 


8b 


180.4 


192.6 


0.5 


6.6 


2.5 


8.7 


9 


181.9 


197.1 


2.0 


11.1 


4.0 


13.2 


10 


177.9 


183.9 


-2.0 


-2.1 


0.0 


0.0 



TABLE IV: Values of the AIC, BIC, AAICia and 2 In Bsi (where s is the index of the reference model) for the flat models from 
Table El 



where flk.o = — 





case 


AIC 


BIC 


AAIC^i 


2 In Bu 


AAIC,2 


2\nB2^ 




1 


178.5 


187.7 


0.0 


0.0 


-2.4 


0.7 




2 


180.9 


187.0 


2.4 


-0.7 


0.0 


0.0 




3 


180.8 


189.9 


2.3 


2.2 


-0.1 


2.9 




4 


327.0 


336.2 


148.5 


148.5 


146.1 


149.2 




5 


182.1 


194.3 


3.6 


6.6 


1.2 


7.3 




6a 


180.3 


192.5 


1.8 


4.8 


-0.6 


5.5 




6b 


183.8 


196.0 


5.3 


8.3 


2.9 


9.0 




7 


180.3 


189.4 


1.8 


1.7 


-0.6 


2.4 




8 


186.6 


195.8 


8.1 


8.1 


5.7 


8.8 




9 


183.9 


196.2 


5.4 


8.5 


3.0 


9.2 




10a 


180.1 


192.3 


1.6 


4.6 


-0.8 


5.3 




10b 


180.2 


192.4 


1.7 


4.7 


-0.7 


5.4 


7?f^and 
















:f{x)-- 


= sinh(a;) 








for 






H^)-- 


= X 








for 






H^)-' 


= sin(a;) 








for 







fc<0 
fc = 
fc>0. 



Table IIIII gives the value of the AIC and BIC for flat models from Table U It also contains values of AAIC^^ and 
2\n Bsi, where model s is our reference model: the ACDM model (indexed by 1) in the first case and model which 
minimize both the AIC and BIC quantities - the flat FRW model with noninteracting the Chaplygin gas and baryonic 
matter (indexed by 10) in the second one. These values are also illustrated in Figures lA, IB and Figures 2A, 2B 
respectively. Table Hvl cont ains the same quantities for flat models from Tabled Here the first reference model is that 
which minimizes the AIC quantity ~ the Cardassian model (indexed by 1), the second one is that which minimizes 
the BIC quantity - the DGB model (indexed by 2). Note that in this case these models are different (on the contrary 
with situation in Table Q] where the same model minimizes both the AIC and BIC indices). Figures 3 A, 3B and 
Figures 4A, 4B illustrate values of AAIC^s and 2\n Bgi respectively. Additionally we compare all flat models from 
Table HTl with the ACDM model (those comparisons are illustrated on Figures 5 A and 5B.) 

In Figures 6A, 6B are presented values of AAICis^;,, and 2lnBs^jci respectively for set of brane models from 
Table Ull {2, 3,4, 5, 6a, 66}. Here the reference model for the AIC analysis is the Shtanov Brand model (indexed by 
6a), which minimizes the AIC quantity in the set under consideration, whereas the DGB model (indexed by 2) for 
the BIC analysis (this one minimizes the BIC quantity in the models set). 
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FIG. 1: Value of A) AAICii = AICi - AICi and B) 2ln Bu for models from Table|I| 
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FIG. 4; Value of A) AAICi2 = AICi - AIC2 and B) 2lnB2i for models from TableEl 
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ence model-the model 10 in both cases. 



TABLE V: Values of the prior and posterior probability for models from Table |I] 



case 



7 



8a 



8b 



9 



10 



prior 
posterior 



0.0909 0.0909 0.0909 0.0909 0.0909 0.0909 0.0909 0.0909 0.0909 0.0909 0.0909 
0.1315 0.0293 0.0240 0.3399 0.0561 0.0308 0.0025 0.0046 0.0048 0.0005 0.3757 



In Figure 7A are illustrated values of AAICis for models which have substantial support with respect to the reference 
model in Figures 2A and Figure 3A: set of models from Table HJ {1,2,4,5,6,10} together with set of models from 
Table im {1, 6a, 7, 10a, 10b}. Here the reference model is this which has a minimal value of the AIC quantity in the set 
of models under consideration - the flat FRW model with noninteracting Chaplygin gas and baryonic matter (model 
indexed by 10 from Table P). In Figure 7B are presented values of 2\n. Bsi for models from Figures 2B and 4B for 
which evidence against them with respect to the reference model is not worth more than a bare mention: Table Q] 
{1, 4, 10} and Table Hll {1, 2}. Here the reference model is this which minimizes the BIC quantity ~ the model indexed 
by 10 from Table U 

In Tables IVI IVll IV 111 and I Villi are gathered values of prior and posterior probabilities (see (|HJl) for models from 
Table U Table ITU set of brane models from Table ITll set of the best models from Table |2 and Table HTl respectively. 

III. CONCLUSION 

Main aim of the paper was exploring the Bayesian framework of model selection into discussion which cosmological 
model describe present accelerating phase expansion of the Universe. In principle there are two types of explanation 
why current Universe is accelerating. In the first group it is postulated existence of perfect fluid which violate 
the strong energy condition. The nature (Lagrangian) of such matter called dark energy is unknown although the 
cosmological constant is the most popular candidate for dark energy description. 

The second group of explanations is based on hypothesis that dark energy could actually be the manifestation 
of a modification to the Friedmann equation arising from 'new physics' (e.g. extra dimensions, generalized general 
relativity, etc.). Calculating the corrections to the standard FRW equation we can explore the phenomenology 
and different evolutional scenarios. Therefore there is no single hypothesis, instead there are several well-supported 
hypotheses (i.e. dark energy models) that are being entertained. It is just realization of concept of 'Multiple Working 



TABLE VI: Values of the prior and posterior probability for models from Table ITTl 



6a 6b 7 



9 10a 10b 



prior 
posterior 



0.083 0.083 0.083 0.083 0.083 0.083 0.083 0.083 0.083 0.083 0.083 0.083 
0.122 0.173 0.041 0.000 0.004 0.011 0.002 0.052 0.285 0.285 0.012 0.012 
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TABLE VII: Values of the prior and posterior probability for set of brane models from Table IITI 



prior 
posterior 



6a 6b 



0.167 0.167 0.167 0.167 0.167 0.167 
0.749 0.175 0.019 0.000 0.047 0.008 



TABLE VIII: Values of the prior and posterior probability for set of models from Table Q {1, 4, 10} and Table ITTI {1, 2}. 



prior 
posterior 



10 



0.20 0.20 0.20 
0.13 0.34 0.38 



0.20 0.20 
0.06 0.08 



Hypotheses' advocated by Chamberlin ^j. In Tables III and Hll we completed the 20 models in two classes which 
cover both types of explanation of acceleration of the Universe. Then we adopt the model selection methods to 
obtain 'strength of evidence' comparison and ranking of candidate hypotheses of dark energy. Providing quantitative 
information to judge the 'strength of evidence' is a key point of ranking models. The hypothesis testing (which 
only provides qualitative information significant vs. nonsignific ant) is particularly limited in the model selection (for 
discussion of the likehhood based strength of evidence see Ref. |53). 

Hence we obtain the set of models which are recommended by the AIC and Bayes factor. The results are following: 

1. While the AIC recommended dark energy models which belong to the set {10,4,5,6, 1,2} the Bayes factor is 
more restrictive because recommended models which are elements of subset {10,4, 1}. This is due to different 
values of the coefficient in the definition of AIC and BIC quantities, so called penalty term for more complex 
models. In the AIC definition this coefficient is always equal 2 whereas in the BIC definition it depends on 
simple size (here it is nearly equal 5). 

2. Analogous recommendation can be performed for the class of models with a modified Friedmann equation, 
namely the AIC recommends a set of models {1, 8, 10a, 106, 6a, 7} and the Bayes factor favors class of models 
{2, 1} which has non-empty intersection with models recommended by the AIC. 

3. One can construct also ranking within the best models in each category I and II by AIC and Bayes factor 
respectively. Then we obtain a set of models recommended by the AIC {10,4, 5, 6, 1, 2}/ U {1}// and by the 
Bayes factor {10, 4, 1}/. 

Note that both quantities pointed out that model with noninteracting Chaplygin gas and baryonic matter is 
the best one from both sets of models under consideration. This is due to more a priori information about this 
model which we include in calculation. Two model parameters were fixed here: ri,„.o smd m. In situation when 
both parameters are fitted the conclusion is changed. All additional information which we include in calculation 
are necessary and can change our final inference. 

It should be point out that dark energy models better explain SNIa data than models with a modified Friedmann 
equation. 

4. We can perform ranking of the models within 'brane paradigm' and then we obtain models {6a, 3, 2, 5} recom- 
mended by the AIC and model {2} preferred by the Bayes factor. 

5. One can conclude that the fiat ACDM model (model indexed by 1) has still substantial support with respect 
to better models. Note that the AIC indicates that there is no difference between the ACDM {k = 0) model 
and the ACDM (A; ^ 0) model (indexed by 2), both of them fit the data equally well. Whereas the Bayes factor 
denotes that there is an evidence against model 2 with respect to model 1. 

Note that Bayes factor indicates that the fiat ACDM model better explains SNIa data than models gathered in 
Table II. 

6. For completeness we calculate posterior probability which measures probability for the model being the best 
one among the class of models under consideration (being the most favored model by data in hand). Then one 
can observe how prior believe about model probability change after inclusion data to analysis (see Tables IVl I VII 
I VIII IVIII|I . Note that for all models recommended by the BIC posterior probability is greater than the prior 
one, whereas for models which do not belong to a preferred set of models the posterior is smaller than the prior. 
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